Symmetry breaking for ratchet transport in presence of interactions and magnetic 

field 
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We study the microwave induced ratchet transport of two-dimensional electrons on an oriented 
semidisk Galton board. The magnetic field symmetries of ratchet transport are analyzed in pres- 
ence of electron-electron interactions. Our results show that a magnetic field asymmetric ratchet 
current can appear due to two contributions, a Hall drift of the rectified current that depends only 
weakly on electron-electron interactions and a breaking of the time reversal symmetry due to the 
combined effects of interactions and magnetic field. In the latter case, the asymmetry between pos- 
itive and negative magnetic fields vanishes in the weak interaction limit. We also discuss the recent 
experimental results on ratchet transport in asymmetric nanostructures. 

PACS numbers: 72.40. +w,73.63.-b,05.45.Ac 
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I. I. INTRODUCTION 

The appearance of a directed flow induced by a zero- 
mean monochromatic force is a generic nonequilibrium 
phenomenon known as the ratchet effect (see e.g. [3-[3])- 
The ratchet transport at nanoscale attracts a significant 
interest in last years (see and Refs. therein). The 
electron ratchet currents induced by a monochromatic 
oc-driving have been experimentally observed in asym- 
metric mesoscopic structures at high Q , very low @ and 
GHz Q frequencies. In the later experiments the ratchet 
effect was visible even at room temperatures but these 
experiments were lacking a detailed analysis of magnetic 
field resonance effects, typical of two-dimensional elec- 
tron gas (2DEG) in a magnetic field @, H| and a depen- 
dence study of ratchet current on microwave polarization. 
We note that the presence of resonances in a resistivity 
dependence on a magnetic field and polarization depen- 
dence of ratchet transport ensure that a voltage induced 
by a microwave radiation appears due to an asymmetry 
of lattice structure and not due to other asymmetries 
potentially present in an experimental devices. 

The theoretical studies of 2DEG deterministic ratchet 
transport on a semidisk Galton board has been started 
in [lcj and further extended in 0, El, 13 0- The 
theoretical studies show that the ratchet effect exists 
not only for 2DEG but also for electrons in graphene 
plane with oriented semidisk lattice [l5j . The nu- 
merical simulations and analytical theory developed in 
[T5L ITl| . fl2l. IT3jj . EEj have been done for noninteract- 
ing electrons. They established that a directed chaotic 
deterministic ratchet transport emerges in such asym- 
metric nanostructures due to microwave radiation. The 
direction of ratchet current can be controlled by the ra- 
diation polarization. The theoretical studies show that 
the ratchet effect also exists in a presence of moderate 
and strong interactions between electrons fl6j . 

The theoretical works initiated the detailed experimen- 



tal studies of the chaotic deterministic ratchet 2DEG 
transport on the semidisk Galton board of antidots per- 
formed by Grenoble group fl~7j . These experiments 
clearly demonstrated the existence of ratchet transport 
in a high mobility 2DEG based on AlGaAs/GaAs hetero- 
junstions with semidisk array. The polarization depen- 
dence of ratchet current is found to be in a qualitative 
agreement with the theory dependence for noninteract- 
ing electrons. It is also experimentally shown [17| that 
the ratchet is absent in arrays with circular antidots en- 
suring that the effect is produced by semidisks and not 
by always present device asymmetries. More recently 
the Grenoble group performed ratchet experiments with 
2DEG in Si/SiGe heterostructures [30 where the in- 
teraction effects between electrons are expected to play 
more important role (l4j |. The characteristic feature of 
these experiments is the dependence of ratchet transport 
on a magnetic field in presence of interactions. Thus it 
is important to understand the properties of ratchet of 
interacting particles in a magnetic field which creates a 
symmetry breaking in space and time. We note that pre- 
viously the theoretical investigations were done only for 
noninteracting electrons in a magnetic field EEEII or f° r 
interacting electrons without magnetic field [16j . Thus in 
this work we perform a more general study analyzing the 
properties of ratchet transport in presence of interactions 
and magnetic field. 

Indeed, an applied magnetic field induces a chiral 
movement of charge in a conducting sample. This chi- 
rality may be revealed in optical measurements such as 
the Faraday effect [1(| EJ. An emergence of a static 
magnetization due to an oc-electric field is known as 
the inverse Faraday effect (see e.g. [11). Hence, one 
can expect that transport properties of a chiral structure 
will strongly depend on the sign of magnetic field. How- 
ever, this argument, based on spatial symmetries, should 
also take into account that the equilibrium transport 
measurements are also constrained by the time rever- 
sal symmetry that implies Onsager-Casimir reciprocity 
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relations [23|, [24J. Due to these relations a two termi- 
nal conductance is always symmetric with magnetic field 
masking the chirality. However, since the time reversal 
symmetry is valid only for equilibrium samples, it can 
be destroyed for measurements in a nonlinear transport 
regime. The appearance of magnetic field asymmetry 
in a nonlinear two terminal transport, has attracted a 
significant theoretical and experimental attention in the 
mesoscopic physics community. Indeed, a microscopic 
disorder potential has no symmetry, and the absence of 
self- averaging in coherent samples makes such investiga- 
tions possible. Surprisingly, even if all symmetries are 
broken in the out of equilibrium regime theoretical cal- 
culations predict that interactions are required to observe 
a magnetic field antisymmetric component in nonlinear 
conductance. While measurements on coherent samples 
seem to support these predictions, measurements on sam- 
ples with artificial asymmetric structures may produce 
asymmetric (or even antisymmetric) nonlinear transport 
even in regimes where interactions do not seem to play 
a role (for e.g. high density) [13, IHl- Thus, the inves- 
tigations of ratchet transport in presence of interactions 
and magnetic field will allow us to analyze the Onsager- 
Casimir reciprocity relations in a new frame. 

In order to gain a deeper insight on the symmetry prop- 
erties of nonlinear transport, we investigate numerically 
an interacting two-dimensional electron gas (2DEG) with 
a periodic array of asymmetric (semidisk) antidots ori- 
ented in a preferential direction u sing the square lattice 
of oriented semidisks discussed in jl5l. Il6|. In this model 
a semidisk of radius is placed in a square of size R 
and then this square covers periodically the whole (x,y) 
plane (see inset of Fig. [T]). We consider a homoge- 
neous monochromatic linearly polarized electric 
field E = E a coswi(cos£?, smff). This field creates a recti- 
fied current flow due to the asymmetric structure of the 
antidot superlattice. The direction of the flow can be 
controlled by the polarization of the microwave field and 
by a magnetic field perpendicular to the 2DEG plane. 

The quantitative description of the ratchet current can 
be obtained on the basis of kinetic theory 14j . However 
this theory rapidly breaks down in presence of a magnetic 
field since it does not capture the particle dynamics when 
cyclotron radius becomes of the order of the antidot ra- 
dius. Moreover it was shown recently [l6| that interac- 
tions could strongly modify the rectified current, hence 
this system allows to investigate in detail the role of inter- 
actions on the magnetic field symmetry properties of the 
nonlinear response. In fact, we show that the antisym- 
metric component of the rectified current can be splitted 
in two terms. One term is weakly dependent on interac- 
tions and can be interpreted as a Hall drift of the rectified 
current in the preferential direction fixed by the semidisk 
superlattice. The second term vanishes in absence of in- 
teractions as predicted by theories 26], [27], [28} . Hence, 
depending on the measurement geometry, the antisym- 
metric component of the rectified current may vanish or 
not in absence of interactions. 



II. II. MODEL DESCRIPTION AND RESULTS 

In this work we consider a 2DEG with elastic semidisk 
scatterers of radius r^, oriented in direction e x = x, and 
placed in a periodic square lattice of size Rx R (see inset 
of Fig [1]) . The electron motion is affected by an electric 
microwave field E cos tot — Eq (cos 9, sin 9,0) cos uit linear 
polarized at angle 9 to e x , and a transverse, uniform and 
constant magnetic field B oc e z /Rr, (inversely propor- 
tional to the Larmor radius Rl)- The system also inter- 
acts with a Nose-Hoover thermostat which equilibrates 
the ensemble of particles to the Boltzmann distribution 
with temperature T = in a characteristic time th 

(see e.g. [29j|). The electron interactions are treated in 
the frame of the mesoscopic multi-particle collision model 
proposed by Kapral (see e.g, H(|). The method consists 
in dividing the coordinate space of each square of size 
Rx R with N particles in N ce i collision cells. Inside each 
cell the collisions are mimic with a rotation of all particle 
velocities on a random angle in the moving center-of-mass 
frame, preserving the total momentum and energy of the 
system. These rotations are done with a repetition pe- 
riod given by a characteristic time tk- In this way large 
and small values of tk correspond to weak and strong 
interactions respectively. In this work we follow our pre- 
vious studies of interactions effects on ratchet transport 
[l6| where the interactions were treated in the frame of 
Kapral approach. 
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FIG. 1: Dimensionless average ratchet velocity v v /vt in di- 
rection y as a function of a magnetic field B given by dimen- 
sionless parameter r^jRh oc B. The electric field is linearly 
polarized in y axis (8 = tt/2) with amplitude and frequency 
-Bo = 0.1 and uj = 0.1, and the temperature of Hoover ther- 
mostat is T = 1 with th = 10. The interactions between 
electrons goes from strong interaction regime {tk /th = 0.02) 
to weak interactions (tk/th = 0.5) as it is shown on legend. 
The inset panel shows the directions of electric and magnetic 
fields in a general case, and the distribution of semidisk peri- 
odic array with an orientation direction in x axis. The geo- 
metric ratio is fixed at R/r c i — 4 as it is shown on the inset. 
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We fix the geometric ratio R/rd = 4 in order to work at 
low antidot density and to avoid geometrical particulari- 
ties which may exist for R/rd ~ 1. The number of parti- 
cles is fixed at N = 10 4 , and for numerical simulations we 
choose dimensionless parameters T — 1, rd — I, th = 10, 
and electron charge and mass e = m e = 1. The Kapral 
grid for interactions is fixed to have N ce i = 100 x 100 
cells in the whole space region R 2 . Therefore the only 
parameter controlling the interaction strength is the in- 
verse Kapral time I/tr. This choice of parameters is 
similar to those used in [l6j]. We remind that tjj deter- 
mines the relaxation time to the equilibrium. 
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FIG. 2: Dimensionless average ratchet velocity in direction 
x (v x /vt) as a function of magnetic field. Top panel shows 
the strong interactions regime (with tk/th = 0.02) for the 
rest of the parameters specified in caption of FigQ] The elec- 
trical field is polarized near x axis, where dashed red curve 
represents the symmetrical case of = 0. The small devi- 
ation angles are shown at 9 — 7r/40 by dotted green curve, 
9 — n/20 by solid blue curve, and 9 — ir/10 by dot-dashed 
black curve. Bottom panel shows the case of fixed polarization 
angle 9 = 7r/20, and different strength of interactions between 
particles. The interaction times are tk /th = 0.02 (solid blue 
curve), tk/th = 0.03 (dashed red curve), tk/th = 0.1 (dot- 
ted green curve) and tk/th = 0.5 (dot-dashed black curve). 
Inserts in both panels show the dimensionless asymmetry 
[v x (r d /RL) — v x (—r d / 'Rl)]/vt as a function of r d /RL for the 
same parameters marked by color. Orange arrows show the 
values of magnetic field analyzed in Fig. [3] 



The steady-state net current of the system is propor- 
tional to the average particle velocities. This current can 
be described as a vector j splitted in two components of 
different magnetic field symmetry: j = j s + j a where j s 
are symmetric and j a antisymmetric components. These 
two vectors depend quadratically on the alternating elec- 
tric field amplitude E. and can also depend on e x and 



B, therefore their most general expression reads: 

j s = / 1 (B)(e a! -E)E + / a (B)E a e a! (1) 
j„ = gi(B)(e x -E)(BAE)+g 2 (B)E 2 (BAe x ) 

Here fi(B), gi{B) (i = 1,2) are functions of the modulus 
of the magnetic field B = |B| depending implicitly on 
the rest of model parameters. Specifically, we focus on 
the dependence of the antisymmetric component of the 
current on the strength of interactions between electrons. 

We note that both terms of second line of Eq. [T]can be 
decoupled for different values of 8. For electric and mag- 
netic fields used in this model, the dimensionless velocity 
can be written as 
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- <7i cos 2 8 + g 2 



where fi and gi (i = 1,2) are symmetric and anti- 
symmetric functions of the magnetic field in z direc- 
tion (or related Larmor radius R^ — vt/B). These 
functions are proportional to the functions in Eqs. Q] 
fi oc fi{B) and gi oc gi(B)/RL (where we take Rl > 
and Rl < for magnetic field in z and — z directions re- 
spectively). Following Eqs. [2] we can note that for 6 = 
and 6 = tt/2 velocities in x direction are symmetric while 
the ^-components are antisymmetric. 

In the case of polarized electric field in y direction 
(8 = 7r/2), v y can be asymmetric in the magnetic field 
only due to the second term of second Eq. [2] noted as 
cj2- This case is illustrated in Fig. [1] for the dimensionless 
quantities v v /vt plotted versus rd/RL at different values 
of Kapral time tk and the rest of parameter specified 
in the caption. The data of Fig. Q] show that the depen- 
dence of ratchet velocity v v /vt is an asymmetric function 
of TdjRh B with maximum of \v y \ at TdjRh ~ ±0.15. 
The amplitude of this maximum increases by a factor 2 
with the increase of interactions (decrease of tk)- 

Following the first line of Eq.[5]the asymmetry given by 
gi can be analyzed via parameter dependence of v x /vt 
at small angles of linear polarization of electric field di- 
rected both along x and y. The emergence of g\ from a 
symmetric behavior in the magnetic field is shown in top 
panel of Fig. [2J for v x /vt with small values of 8 in the 
strong interactions regime (tk/th = 0.02). In the case 
of 8 = it/ 20, bottom panel of Fig. [2] shows the behavior of 
asymmetry in v x /vt for different interaction times, going 
from tk/th = 0.02 to tk /th = 0.5. The striking feature 
of bottom panel of Fig. 2 is that the asymmetry in a mag- 
netic field is rather strong at strong interactions while in 
the limit of weak interactions the asymmetry completely 
disappears and we recover the symmetric curve as a func- 
tion of a magnetic field. 

A deeper analysis of the symmetry properties can be 
done by considering the flux of velocities in coordinate 
space. It is known that for interacting particles the av- 
erage velocity behavior can be rather complex with an 
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emergence of some vortexes [16]. In Fig. [3] we present 
the velocity flux for weak and strong interactions regimes 
with positive and negative magnetic fields. The analyzed 
values of Larmor radius correspond to the relative min- 
ima of v x /vt marked with arrows in Fig. [2] at negative 
and positive values of rd/Rh- Fig.[3]shows the flow struc- 
ture at Tdj ' Rl — 0.16 on top panels and at rd/Rh — —0.1 
on bottom panels. The polarization angle in four panels 
is fixed at = 7r/20, and therefore the reflection symme- 
try y — » — y is not preserved. The data show that the 
vortex asymmetric structure is more pronounced in the 
case of strong interactions on right panels. 





FIG. 4: (color online) Polarization dependence of x- 
component of the averaged ratchet velocity v x /vt, for weak 
interactions between particles in top panel (tk/th = 0.5), 
and for strong interactions in bottom panel (jk/th = 0.02). 
Four different values of magnetic fields are shown in both 
panels: zero magnetic field (rd/Rh = oo, solid blue curve), 
rd/Rh = 0.053 by dotted red curve, rd/Rh = 0.133 by dashed 
green curve and rd/Rh = 0.53 by dot-dashed black curve. 



the data of Fig.|4]also show that the asymmetry of ratchet 
transport appears only in a presence of interactions. 



III. III. DISCUSSION 



FIG. 3: Map of local averaged velocities in (x, y) with x, y £ 
[—3,3) and linear polarized electrical field 6 = 7r/20 (shown 
by green line inside semidisk scatterer); other parameters are 
the same as in Fig. Q] Fhe left panels show the cases of weak 
interactions between particles at tk/th = 0.5, while the right 
panels show strong interactions regime at tk/tb = 0.02. The 
values of magnetic field are given by rd/Rh = 0.1591 on top 
panels, and rd/Rh = —0.0955 on bottom panels (pointed with 
orange arrows in Fig. [2]). The velocities are shown by arrows 
which size is proportional to velocity amplitudes, which is also 
indicated by color [from yellow (gray) for large to blue (black) 
for small amplitudes]. 

The polarization dependence of ratchet current in 
a;— direction is analyzed in Fig. @] We see that even for 
various magnetic fields the dependence on polarization 
angle 9 is essentially symmetric at weak interactions (top 
panel) while at strong interactions (bottom panel) we 
have a strongly asymmetric behavior at moderate mag- 
netic fields {rd/Rh = 0.053,0.133). At a relatively large 
magnetic field (rd/Rh — 0.53) the amplitude of ratchet 
current becomes rather small since the Larmor radius be- 
comes smaller than the distance between semidisks. Thus 



In this work we study the symmetry properties of the 
ratchet transport on an oriented semidisk antidot super- 
lattice. We show that while the ratchet flow on the 
oriented semidisk superlattice (along a;— axis) does not 
depend on the sign of magnetic field for a noninteract- 
ing 2DEG, interactions can give rise to an antisymmet- 
ric component of the flow in the semidisk direction as a 
function of a magnetic field. This result is consistent 
with the case of mesoscopic samples where deviations 
from Onsager-Casimir reciprocity relations were shown 
to occur only in presence of electron-electron interactions 
26], 27], 28]. On the contrary, the flow perpendicular 
to the semidisk direction (y— axis) is asymmetric even 
when interactions are absent. We argue that the origin 
of this component of the flow is a Hall drift of the rec- 
tified current. This contribution is absent in disordered 
mesoscopic samples because the translational symmetry 
is broken by a disorder potential where no preferential 
direction is present for the flow. 

The comparison with [13], 18|,[IU highlights various 
aspects of these skillful experiments. At weak electron- 
electron interactions typical of AlGaAs/GaAs herero- 
junctions [17] there is a qualitative agreement between 
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the theory and experiment on polarization dependence 
and approximately symmetric current response on sign 
change of a magnetic field, even if the absolute quanti- 
tative values differ from theoretical predictions (see [l6[ 
for a more detailed discussion) . For experiments with 
Si/SiGe heterostructures [HI, |T^| the effects of interac- 
tions are stronger and asymmetric response in a magnetic 
field is clearly observed in experiments (see e.g. Fig. 3 in 
[l8|). In this case the change of polarization from 9 = 
to 9 = ir/2 does not change the sign of photovoltage that 
also appears in theoretical models at strong interactions 
as discussed in [h|- Also the ratchet effect is clearly 
suppressed in experiment and theory at large magnetic 
fields. However, at the same time there are significant 
differences between experiments [lH [l9| and theoretical 
results presented here. Indeed, for 9 = ir/2 the theory 
predicts a change of the photovoltage sign upon inversion 
of the magnetic field (see Fig. Q] above) while there is no 
such sign change in experiments (see e.g. Fig. 3 in (l8j 
and in [19( ) ■ It is possible that finite sample size leads to 
a certain charge accumulation in the experimental setup 
(see indications for that in Fig. 4 in 0]) that may explain 



the difference with the theory where the analysis is done 
for an infinite lattice size. We should also note that the 
present studies were done for a square lattice of semidisks 
while in the experiments [17H19J were performed on an 
hexagonal lattice of semidisks. However, in both cases 
the density of semidisks is relatively low (since R/r<i = 4 
here and R/rd ~ 5 in 0, 0]). Thus, the ratchet trans- 
port is created mainly by scattering on a one semidisk 
(see discussion in 0) and the difference between square 
and hexagonal lattices is not expected to be important. 

The discussed caparison between the present theoreti- 
cal studies and the most advanced experiments reported 
in [13)050 shows that the further experimental and 
theoretical research of the electron ratchet transport in 
asymmetric nanostructures with dynamical chaos repre- 
sents a significant fundamental scientific interest. 

We thank M.L. Polianski for fruitful discussions on out 
of equilibrium transport in mesoscopic systems. This 
work was supported in part by the French PNANO ANR 
project NANOTERRA, one of us (A.D.C.) acknowledges 
support from St Catharine's College. 



[1] J. Jiilicher, A. Ajdari and J. Prost, Rev. Mod. Phys. 69, 
1269 (1997) 

[2] R.D. Astumian and P. Hanggi, Phys. Today 55(11), 33 
(2002) 

[3] P. Reimann, Phys. Rep. 361, 57 (2002) 
[4] P. Hanggi and F. Marchesoni, Rev. Mod. Phys. 81, 387 
(2009) 

[5] A. Lorke, S. Wimmer, B. Jager, J. P. Kotthaus, 
W. Wegscheider, and M. Bichler, Physica B 249-251, 
312 (1998) 

[6] H. Linke, T.E. Humphrey, A. Lofgren, A.O. Sushkov, 
R. Newbury, R.P. Taylor, and P. Omling, Science 286, 
2314 (1999) 

[7] A.M. Song, P. Omling, L. Samuelson, W. Seifert, I. Sho- 
rubalko, and H. Zirath, Appl. Phys. Lett. 79, 1357 (2001) 

[8] D. Weiss, M. L. Roukes, A. Menschig, P. Grambow, K. 
von Klitzing, and G. Weimann, Phys. Rev. Lett. 66, 2790 
(1991) 

[9] GM. Gusev, Z.D. Kvon, V.M. Kudryashov, L.V. Litvin, 
Y.V. Nastaushev, V.T. Dolgopolov, and A. A. Shashkin, 
JETP Lett. 54, 364 (1991) [Pis'ma Zh. Eksp. Teor. Fiz. 
54, 369 (1991)]. 

[10] A. D. Chepelianskii, and D. L. Shepelyansky, Phys. Rev. 
B 71, 052508 (2005) 

[11] G. Cristadoro, and D.L. Shepelyansky, Phys. Rev. E 71, 
036111 (2005) 

[12] A. D. Chepelianskii, Eur. Phys. J. B, 52, 389 (2006) 
[13] M.V. Entin, and L.I. Magarill, Phys. Rev. B 73, 205206 
(2006) 

[14] A.D. Chepelianskii, M.V.Entin, L.I.Magarill and 
D.L. Shepelyansky, Eur. Phys. J. B 56, 323 (2007) 

[15] L.Ermann and D.L. Shepelyansky, Eur. Phys. J. B 79, 
357 (2011) 



[16] A.D. Chepelianskii, M.V.Entin, L.I.Magarill and 

D.L. Shepelyansky, Phys. Rev. E 78, 041127 (2008) 
[17] S. Sassine, Yu. Krupko, J.-C. Portal, Z.D. Kvon, R. Mu- 

rali, K.P. Martin, G. Hill and A.D. Wieck, Phys. Rev. B 

78, 045431 (2008) 
[18] I. Bisotto, E.S. Kannan, R. Murali, T.J. Beck, L. Jalabert 

and J.-C. Portal, Nanotechnology 22, 245401 (2011) 
[19] E.S. Kannan, I. Bisotto, J.-C. Portal, R. Murali and 

T.J. Beck, Appl. Phys. Lett. 98, 193505 (2011) 
[20] M. Faraday, Faraday's Diary, vol.IV, Nov 12, 1839 - 

June 26, 1947, Th. Martin Ed., George Bell and Sons, 

London (1933). 

[21] L.D. Landau and E.M. Lifshitz, Electrodynamics of con- 
tinuous media, Pergamon Press, New York (1960) 
[22] M.L. Polianski, Phys. Rev. B 80, 241301 (R) (2009) 
[23] L. Onsager, Phys. Rev. 37, 405 (1931); ibid. 38, 2265 
(1931) 

[24] H.B.G. Casimir, Rev. Mod. Phys. 17, 343 (1945) 

[25] A.D. Chepelianskii and H. Bouchiat, Phys. Rev. Lett. 

102, 086810 (2009) 
[26] D. Sanchez and M. Biittiker, Phys. Rev. Lett. 93, 106802 

(2004) 

[27] B. Spivak and A. Zyuzin, Phys. Rev. Lett. 93, 226801 
(2004) 

[28] M.L. Polianski and M. Biittiker, Phys. Rev. Lett. 96, 
156804 (2006) 

[29] W. G. Hoover, Time Reversibility, Computer Simulation, 

and Chaos, World Sci., Singapore (1999) 
[30] A. Malevanets and R. Kapral, in Novel Methods in Soft 

Matter Simulations, Eds. M. Karttunen, I. Vattulainen, 

and A. Lukkarinen, Lecture Notes in Phys. (Springer, 

Berlin) 640, 116 (2004) 



